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Abstract 

' We analyze a finite element/boundary element procedure to solve a non- 

^ | . convex contact problem for the double-well potential. After relaxing the 

associated functional, the degenerate minimization problem is reduced to 
a boundary/domain variational inequality, a discretized saddle point for- 
mulation of which may then be solved numerically. The convergence of 
the Galerkin approximations to certain macroscopic quantities and a cor- 
responding a posteriori estimate for the approximation error are discussed. 
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m ■ 1 Introduction 

CN ' Adaptive finite element / boundary element procedures provide an efficient 

and extensively investigated tool for the numerical solution of uniformly elliptic 
transmission or contact problems. However, models of strongly nonlinear mate- 
rials often lead to nonelliptic partial differential equations, where the standard 
Hilbert space techniques are no longer appropriate to analyze the computational 
methods. In a previous work [8] we showed that certain mixed I? — L p -Sobolev 
spaces provide a convenient setting to study contact problems for monotone 
operators like the p-Laplacian. This article extends the approach to nonconvex 
functionals, discussing the prototypical model case of a double-well potential 
in Signiorini and transmission contact with the linear Laplace equation. As a 
proof of principle, it intends to clarify the mathematical basis - including well- 
posedness, convergence, a priori and a simple a posteriori estimate - of adaptive 
finite element / boundary element methods in this highly degenerate nonlin- 
ear setting. The methods readily extend to certain systems of equations from 
nonlinear elasticity, frictional contact or more elaborate a posteriori estimates. 

Let SI C M. n be a bounded Lipschitz domain and d£l = Tt U T s a decom- 
position of its boundary into disjoint open subsets, T t ^ 0. We consider the 
problem of minimizing the functional 

*(«i,«2)= / W(Vui) + - [ |Vu 2 | 2 - / /ui-(to,«3|«i> 
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with nonconvex energy density W(F) = \F — Fi\ 2 \F — F 2 \ 2 (iq ^ F 2 £l") 
over the closed convex set 

{( Ul ,u 2 ) G W l >\tt) xW^(Q c ) : (ni-n 2 )|r t = u , (u 1 -u 2 )\r s < «o, ^2 G A}, 

£ 2 = {{« G Wf ' c 2 (tt c ) : A W = in ^" 1 ' 2 (^), v = { ^ |2 _ 

The data / G L 4/3 (^), id G W^i 2 (dft) and u G wi 2 (dft) are taken from the 
appropriate spaces. 

Classical exact minimizers of $ satisfy the equations 

-div LW(Vm) = / in 0, Au 2 =0 in fl c , 
u ■ DW(Vui) — d u u 2 = to on dfi, u\ — u 2 =uq on 1^, 
Ui — U2 < uo, ^ • DW(Vu\) <0, 1/ ■ DW(Vu\)(ui — u 2 — uq) =0 on r s , 

+ radiation condition for u 2 at 00. 

Therefore, the minimization problem for $ is a variational formulation of a con- 
tact problem between the double-well potential W and the Laplace equation, 
with transmission (Tj) and Signiorini (r s ) contact at the interface. 

Nonconvex minimization problems of this type arise naturally when a ma- 
terial in Q passes the critical point of a phase transition into a finely textured 
mixture of locally energetically equivalent configurations of lower symmetry, 
the so-called microstructure. Lacking convexity in f2, the minimum of is 
usually not attained. Nevertheless, it is possible and of practical interest to 
extract average physical properties of the sequences minimizing <3?. Examples 
of such quantities include the displacement in the exterior, stresses, the region, 
where minimizing sequences develop microstructure, or also the gradient of the 
displacement away from the microstructure. Crucially for the use of boundary 
elements, the exterior boundary value on the interface is not affected by the 
presence of microstructure. 

The increasingly fine length scale of the microstructure often prevents the di- 
rect numerical minimization, and starting with works of Carstensen and Plechac 
[3j H] computational approaches based on relaxed formulations have been con- 
sidered. Relaxation amounts to replacing the nonconvex functional by its quasi- 
convex envelope, in our setting the degenerate functional 

$**(ui,u 2 )= I W**{Vu x ) + \ I |Vu 2 | 2 - I M-(i ,n 2 |an). 
Jn * Jn c Jq 

If A = |(F 2 — Fi) and B = |(iq + -^2), the convex integrand W** is given by 
the formula (cf. [3]) 

W**(F) = (max{0, \F - B\ 2 - \A\ 2 }) 2 + 4|y4| 2 |F - B\ 2 - A(A(F - B)) 2 . 

The theory of relaxation for nonconvex integrands shows that the weak limit 
of any ^-minimizing sequence minimizes <!>**. Macroscopic quantities like the 
stress DW** on f2 defined by this weak limit coincide with the averages such as 

2 




the average stress J DW(u) du(u) defined by the Young measure v associated 
to the minimizing sequence. To extract the average physical properties of se- 
quences minimizing <£, it is thus sufficient to understand the minimizers of the 
degenerately convex functional <&**. 

We are thus going to analyze a finite element / boundary element scheme 
which numerically minimizes $** and thereby approximates certain macroscopic 
quantities independent of the particular minimizer. Our approach is based on 
previous works by Carstensen / Plechac O U] and Bartels [I] for double-well 
potentials with Dirichlet or Neumann boundary conditions. Section [2] discusses 
the relaxed problem and identifies several quantities shared by its minimizers. 
A priori error estimates for their computation and convergence are established 
in Section [3) Section H] contains an a posteriori estimate of residual type, on 
which an adaptive grid refinement strategy may be based. 

For later reference, we recall from [3j the following estimates for the relaxed 
double-well potential (E,F € M. n ): 

max{Ci|F| 4 -C 2 ,0} < W**(F) < C 3 + C 4 \F\\ (1) 

\DW**(F)\ <C 5 (1 + |F| 3 ), (2) 

\DW**(F) - DW**(E)\ 2 < C 6 (l + \F\ 2 + \E\ 2 )(DW**(F) - DW**{E))(F - E), 

(3) 

8\A\ 2 \FF - ¥E\ 2 + 2 Qi - F) +^ E) \A(F - E)\ 2 + 2(Q(F) - Q(E)) 2 

<{DW**(F)-DW**{E))(F-E) (4) 

where Q(F) = max{0, \F — B\ 2 — |^4| 2 } and P is the orthogonal projection onto 
the subspace of vectors orthogonal to A. 

2 Analysis of the Relaxed Problem 

We first outline how the minimization problem for $** can be reduced to a 
boundary-domain variational inequality. As it involves the exterior problem, it 
is not affected by the nonconvex part of the functional. See e.g. [8] for a more 
detailed exposition. 

Recall the Steklov-Poincare operator 

S : W^' 2 {dn) -> w~^ 2 {dn), 

a positive and selfadjoint operator (pseudodifferential of order 1, if dQ is smooth) 
with defining property 

d v u 2 \dn = -S(u 2 \an) 
for solutions u 2 € C 2 of the Laplace equation on £l c . Let 

W^ 2 (T S ) = {ve W^' 2 {dQ) : supp v C f s }, X = W 1A (n) x W?' 2 (T S ). 
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Using S, an affine change of variables, 

(ui,u 2 ) h4 (u,v) = (ui - c,u + u 2 \an - ui\an) G X, 
for a suitable c£l reduces the exterior part of $** to T s : 

«2)= / w^*(Vu) + i(5Han+«),u|fln+«)-A(u,«)+C7= J(«,v) + C, 

where 

A(n,u) = (i + 5 , -u ,n| an + w) + / /it 

in 

and C = C(uo,to) is a constant independent of u,v. Therefore, instead of <£** 
one may equivalently minimize J over 

A = {(u, v) G X : v > and (5(tt|fln + « - «o), 1) = if n = 2}. 

A reformulation as a variational inequality reads as follows: Find (u, v) £ A 
such that 

/ DW*%Vu)V(u-u) + (S(u\ dn + v),(u-u)\ 9n +v-v) > \(u-u,v-v) (5) 
for all (u, v) G A 

Convexity and the closedness of A assure that the relaxed functional J 
assumes its minimum. Due to the lack of coercivity, the minimizer may fail to 
be unique, though certain macroscopic quantities are uniquely determined. 

Lemma 2.1. The set of minimizers is nonempty and bounded in X . The stress 
DW**(u), the projected gradient PVu, the region of micro structure {x G £1 : 
Q(Vu) = 0} and the boundary value u\qq + v are independent of the minimizer 
(u,v) G A of J (up to sets of measure 0). 

For the proof, we recall the variant 

INIwM(fi) < ||Vu||i4 (n j + \\u\\ w i i2{rt) (6) 
of Friedrichs' inequality from [8]. 

Proof of Lemma \2.1\ By (pQ) and the coercivity of S, we have 

\J(u,v)\ > C 1 \\V^ HQ) -C 2 vo\n + -Cs\\Mr s +v\\ 2 wh2 ^ ) 

+ 2 Cs ^wh 2 (r t ) ~ ll/lli 4 / 3 (n)ll«llw ,1 - 4 (n) 

— 1 1 tn + S Un 1 1 i \\u\\ l o, , 

11 u U|l iy-2' 2 (9c) n ll w?' 2 (r t ) 



for any minimizer (u, v) € A of J. Consequently 

ll^llf 4^01 + ll^lr. + ^|| 2 i, + ll^ll 2 1, - C||^llivi.4rn , i 
ii ui (»j ii lis V3' a (r s ) Vj' 2 (r t ) ^ J 
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is bounded for some C > 0. The inequality ([6]) easily yields the boundedness 
of\\(u,v)\\x. 

If Vx), (u 2 , 2 ) € A are two minimizers, J is constant on {(ui,Vi) + 
s(u2 — u%,V2 — Vi) : s € [0,1]}: If not, the restriction of J to this set would 
have a maximum > J(ui,v\) = J(?i 2 ,w 2 ) f° r some < s < 1, contradicting the 
convexity of J. Therefore 

(J'(u2,h) ~ J'(ui,Vl), (U2 ~ Ui,V 2 ~ Vl)) = 0, 

or for our particular J 

= [ (DW**(Vu 2 ) - DW**(Vui))V(u 2 -ui) 
Ju 

+ {S((u 2 - ui)\dn + v 2 - vi), (u2 - mi) |an + h - vi)- 

Both of the terms on the right hand side are non-negative, because S is coercive 
and W** convex, and hence 

ui\an + h = u 2 \ 9 Q + h and (DW**(Vu 2 ) - L>iy**(Vui))V(n 2 - «i) = 
almost everywhere. The inequality ([3]), 

\DW**{Vu 2 )-DW**{Vui)\ 2 < (l+\Vu 2 \ 2 MVu 1 \ 2 )(DW*%Vu2)-DW*\Vu 1 ))V(u 2 -ui), 

implies DW**(\7ui) = DW**(\7u2) almost everywhere. The assertions about 
the projected gradients and the region of microstructure are immediate conse- 
quences of inequality flU), 

|PVn 2 - PVni| 2 + (Q(W 2 ) - <9(Vni)) 2 < (DW**(Vu 2 ) - DW**(Vui))V(u2 - 

□ 

In particular, the displacement U2 on $l c is uniquely determined and may 
be computed from u\qq + v with the help of layer potentials. Due to the lack 
of convexity of W, neither u nor Vu needs to be unique. However, Lemma 12. II 
allows to identify subsets of SI, on which these quantities are well-defined. 

Corollary 2.2. a) Let £lt,A be the set of points x G S7 for which the component 
of a hyperplane perpendicular to A through x intersects Y t . Then the displace- 
ment u\q a is independent of the minimizer. 

b) The same holds for the gradient outside the region of microstructure. 

Proof. The proof closely follows the arguments of [3J, Theorem 3. 

a) Let (ui,#i), (u 2 ,v 2 ) G A be two minimizers, and consider w = 112 — u\. 
Because PVui = PVu 2 ) V«) is parallel to A almost everywhere. It is easy to 
see that, therefore, w may be modified on a set of measure zero to yield an 
absolutely continuous function which is locally constant along the hyperplanes 
perpendicular to A. With w\r t being by Lemma 12. 11 w also has to vanish on 
almost every hyperplane hitting Tj. 

b) is a consequence of PVui = PVu 2 , DW**(Vui) = DW**{Vu 2 ) and (@J): 
(Q(Vn 2 ) + Q{Vui))\AV{u 2 - ni)| 2 < (DW**(Vu 2 ) - DW**(X7 Ul ))V(u 2 - ui). 

□ 
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3 Discretization and A Priori Estimates 



We are now going to analyze which quantities can be computed numerically 
with a Galerkin method. 

Let {Th}hei a regular triangulation of Q, C M 2 into disjoint open regular 
triangles K, so that f2 = [j Ke j- h K. Each element has at most one edge on 
(9f2, and the closures of any two of them share at most a single vertex or edge. 
Let hx denote the diameter of K G Th and px the diameter of the largest 
inscribed ball. We assume that 1 < maxxeT h 7^ < R independent of h and 
that h = maxK£T h hx- &h is g° m g to be the set of all edges of the triangles in 
Th, D the set of nodes. Associated to Th is the space 
functions whose restrictions to any K G Th are linear. 

i 2 

d£l is triangulated by {I G Eh '■ I C <9f2}. W^ 2 ' (9f2) denotes the correspond- 

1 2 

ing space of piecewise linear functions, and W h 2 ' (T s ) the subspace of those 
supported on F s . Finally, W~^{d9) C W-^' 2 (dQ), 

A h = An(w h 1 '\n)xwl' 2 (dn)) 

a ndX* = W^(n)xwt 2 (T s ). 



We denote by i h : WjJ' 4 (fi) W 1 ' 4 ^), j h : W^ 2 ' (r s ) W2' 2 (T S ) and 

-- 2 l 

kh '■ W h 2 ' (<9$7) ^ W~z' 2 (d£l) the canonical inclusion maps. A discretization 
of the Steklov-Poincare operator is defined as 

S h = l -{W + (I- K^khiklVkh^kUl - K)) 

from the single resp. double layer potentials V and K and the hypersingular 
integral operator W of the exterior problem. Sh is well-known to be uniformly 
coercive for small h in the sense that there exists h$ > and an /i-independent 
as > such that for all < h < ho 



(S h u h ,u h ) > a s \\u h \ 
Furthermore, in this case 



" S Hl w - i2(aQ) < C S dist w _ h2{dn) (V-\l - K)u, W-^(OCl)) (7) 

for all u G W^> 2 (d£l) and all < h< h . 

As before, (u,v) denotes a minimizer of J over A, while (uh,Vh) minimizes 
the approximate functional 

Jh(u h ,v h )= I W**(Vu h ) + ^{S h (u h \ dn + v h ),u h \an + v h ) - \ h (u h , v h ), 
Jn z 

K{uh,v h ) = (to + S h u ,u h \dn + v h ) + / fu h , 

Jn 

over Ah- For simplicity, abbreviate the stress DW**(Vu) by a and the indicator 
Q(Vu) for microstructure by £. Similarly, write ah and for the corresponding 
quantities associated to %. The following a priori estimate holds. 
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Theorem 3.1. The Galerkin approximations of the stress a, exterior boundary 
values u\qq + v and the other quantities in Lemma \2.1\ converge for h — > 0. 
a) There is an h-independent C > such that 



k-^ll^ (n) + ll(^-^)lan + *-^|| 2 ^, 2(m) 



+ ||PVu - ¥Vu h \\ 2 L2m + + Hh) 1/2 AV(u - fi h )ll£a ( n) + IK " &ll£a 



i 



+ dist^i^^l - K)(n + - no), W h (df>)) 2 . 
For pure transmission conditions, I\ = dSl, the slightly better estimate 



+ii(e + th) 1/2 Av(u - u h )\\ 2 L2(n) + ne - allien) 

C inf {||Vw-V^||?4roi + ||«-^|| 2 i, 



+ dfet w _ ii9 (^(l - *0(fi - no), W h i2 (aO)) 2 



holds. 



Proof. We integrate (J3j) and use Holder's inequality as well as the uniform bound 
on the norm of minimizers (the first assertion in Lemma l2.1j) to obtain 

ll " -0 "^! 2 -!^ ~ / ( ff - - (8) 

Most of the remaining terms on the left hand side are similarly bounded with 
the help of @: 



\¥Vu - PV^||| 2(n) + ||(e + i h ) l/2 AV{u - u h )\\l* m + 11^ - £,h\\h 



(fi) 



< 



((r-<r A )V(«-ti h ). (9) 



Adding the inequalities, 



LHS 2 := Ik-^i^ + IKfi-^lm + C-fiftll^aoj 



+l|pvn - pv« fc ni a(n) + ii(e + ^) 1/2 ^v(^ - %)||| 2(Q) + ne - allien) 

~ /n( a ~ a h)(^u ~ V%) + - u h )\ dn + v -v h ,(u- u h )\dn + v-v h ) 
aVu h - {S{u\ d u + v),u h \ d n + v h ) 



\ a h Vu - {S(u h \dQ + v h ),u\ d Q + v) 

+ / crVn+ (S(u\gn + v),u\on + v) 
Jn 



+ / <ThVuh + (Sh(uh\an + Vh),Uh\dn + Vh) 
+ ((S - S h )(u h \ dn + v h ),u h \ dn + v h ). 
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Let (U, V) £ A, (Uh,Vh) € Ah- Applying the variational inequality ([5]) to the 
third and fourth line and rearranging terms leads to 

LHS 2 < [ aV(U-u h ) + (S(u\ dQ + v),U\ dn + V-u h \ dQ -v h ) 
Jn 

+ / cr h V(U h -u) + (S(u h \ d n + Vh),U h \ dn + V h -u\ dn -v) 
Jn 

+\(u — U,v — V) + \{u h - U h , v h - V h ) 

+ {{S - S h )(uh\da + Vh~ uo), {uh - U h )\ 9 n + v h - V h ) 

aV(U - u h ) + (S(u\aa + v), (U - u h )\ dn + V - v h ) 



+ / aV{U h -u) + (S(u\ dQ + v), (U h - u)\ dQ + V h -v) 



n 



+ / (crh ~ o-)V(C/ /l -u) + {S((u h ~ u)\ 9n + v h - v), (U h - u)\ dn + V h -v) 
Jn 

+\(u — U,v — V) + \{u h - U h , v h ~ V h ) 

+ ({S - S h )(u\ d n + v - u ), (u h - U h )\ dn + v h - V h ) 

+ {{S - S h )((uh - u)\dn +v h - v), (u h - U h )\dn + v h - V h ). 

Holder's inequality tells us that 

{a h - <?)V{U h - u) < \\<Th ~ <r\\ L § {Q) W v ( u h ~ u)\\ L 4(n), 
and the continuity of S allows to bound 

(S((u h - u)\dn + v h - v), (U h - u)\an + V h -v) 
by a multiple of 

\\(u h - u)\ dn + v h - t\\ w ^ m \\(U h - fi)|en +V h - v\\ w ^ {m) 
< e\\(u h - n)\an + v h - 0||^, a(fln) + - £ W h - u)\an + V h - fi||^^ (fln) 
for small e > 0. Similarly the last two lines are, up to prefactors, bounded by 

e\\(u h - u)\an +v h - v\\^ >2{m) + (1 + - £ )\\(U h - &)\ m + V h - 
+ \\(S-Sh)(u\9n + v-n )\\ 2 wh2{m) . 

Thus, choosing (U, V) = (u h ,v h ), 

LHS 2 < hh-a\\ Ll{n) \\V(Uh-u)\\ LH n) 

+Wh ~ Vim + V h - vf wh2{gu) + - S h )(u\ m + v - u )f wl2m 



+ / aV{Uh-u) + (S{u\ dn + v),{Uh-u)\ dn + Vh-v}-\{U h -u,V h -v). 



n 
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If I\ = <9f2, the variational inequality ([5]) becomes an equality, the last line 
vanishes and b) follows. In the general case, we estimate the last line by 

Hl x § (Q) llV(^ - u)\\ Lm + ||/H L 4 (n) ||£4 - fi|| £ 4 (n) 
+ \\S(u\an + v-u )- t \\ w .i, 2(dQ) \\(U h - u)\ m +V h - v\\ w ^ m> 

recalling that 

X(U h -u,V h -v) = {t + Suq, (U h - fi)|an +V h -v)+ [ f(U h - u). 



□ 



In particular, we can stably compute the approximate solutions in the ex- 
terior domain from Uh\dn + Vh- 



4 Adaptive Grid Refinement 

In order to set up an adaptive algorithm, we now establish an a posteriori 
estimate of residual type. It allows to localize the approximation error and leads 
to an adaptive mesh refinement strategy. A related and somewhat more involved 
estimate for the linear Laplace operator with unilateral Signiorini contact has 
been considered in [9]. 

Let (u,v) G A, (ufih) € Ah solutions of the continuous resp. discretized 
variational inequality. We define a simple approximation (iThii, iThv) £ Ah of 
(u, v) as follows: n^u is going to be the Clement interpolant of u, and iThV = Vh- 

The next Lemma collects the crucial properties of Clement interpolation 
(see e.g. 0). 

Lemma 4.1. Let K € Th and E € Eh- Then with ujk = Ux'nx^©-^' an< ^ 
= UwnE^® E ' we have: 

\\u - nhu\\L*(K) ~ hxWuWw^iuK) i 



\u - nh,u\\ L 2(E) 



5. hA \\u\ 



We are going to prove the following a posteriori estimate: 
Theorem 4.2. 

Ik - a h\\j {n) + II (w - Wfc)|an + v- Vhf w i, 2[gn) + ||PVfl - ¥Vu h \\ 2 L 2 {n) 



(n) us. ~ sftiii2(Q) 

< ^ + ^ + r?s + dist w _i, 2(M) (F- 1 (l-K)(^| 9 c + %-n ),^5' 2 (^)) 2 , 
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where 



Vn = J2h K \\f\\ L 4/ 3{K) + h M\["E -<Th]\\LHE) , 

K EndQ=i 

VC = r]c,\ + r]c,2 = ^2 \\{v E ■ a h ) + \\ L 2 {E) + ^ / {y E • <7 ft )_ v h , 
Vs = ^ h 1 ^ 2 \\S h (u h \en + Vh~ u ) + {yaci ■ a h ) - t \\ L 2 (E) . 

EcdQ 

Remark 4.3. a) Also the constant prefactors, suppressed in our notation <, are 
explicitly known. 

b) The main point of this estimate is to show that the a posteriori estimates 
for the contact part ([9]) and the double-well term ([3]) are compatible. More 
sophisticated bounds related to a different choice of 7Th generalize to our setting 
in a similar way. As a simple case, a more considerate (sign-preserving) choice 
of TThV with \\v — ^h^\\L 2 (T s ) ^5 ^ a ||^llL 2 (r a ) could be used to gain an h a in 770,1 
at the expense of modifying 



VC,2 = 




as long as we only assure that f E (iThVh — v) < J E ^\vh for some auxiliary 
interpolation operator tt^ (see e.g. [9J). 

c) As in [8], it is straight forward to introduce an additional variable on the 
boundary to obtain estimates that do not involve the incomputable difference 
Sh — S. Similarly, we might also use the formulation of Bartels pQ with explicit 
Young measures in the interior part. 

Proof of Theorem \4-2\ As in the proof of Theorem 13.14 we start with the in- 
equality 

LHS 2 :=\\a- a h \\\^ + \\(u- u h )\ dn + v - 411^,^ 

+ ||PVn - PVfiftU^n) + m + th) 1/2 AV{u - fi h )||^ (n) + ||e - Zh\\b(n) 
(a - a h )V(u - u h ) + (S((u - u h )\an + v- v h ), (u - u h )\dn + v- v h ). 
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Using the variational inequality and its discretized variant results in 

LHS 2 < X(ii-u h ,v -v h ) - a h V(u-u h ) - (S(u h \ 9n + v h ,(u-u h )\ dn + v -v h ) 

Jn 

= X(u-u h ,v-v h )- a h V(u-u h ) - (S h (u h \m + v h , (u - u h )\ m + v - v h ) 

Jn 

+((S h - S)(u h \ dn + v h , (u - u h )\an + v-v h ) 

< X(u - u h , v - v h ) - I a h V(u-u h ) - (S h (u h \ m + v h , (u - u h )\m + v - v h ) 

Jn 

+ ((S h - S)(u h \an + v h ~ u ), (u - u h )\dn + v - v h ) 

= f(u-U h )- ^ / [ V E ■ Oh]{u - Uh) 

- {Sh{uh\dn + v h - n ) + (van • 07») - to), (u - u h )\an + v - v h ) 



+ I {van ■ 0h) (v ~ v h) 
Jr s 

+((S h - S)(u h \ dn + v h - u ), (u - u h )\ dn + v - v h ) 

for all (uh,Vh) G Ah- Here, ve and fan denote the outward-pointing unit 
normal vector to an edge E C K, resp. to dd, and [ve ■ 07J is the jump of 
the discretized normal stress across E. According to estimate ((7|) for Sh — S 
and Young's inequality, the last term contributes the not explicitly computable 

dis ViW y-1(1 " K) ^ m + ^' w h h2 (9n)) 2 . 

We are going to choose (uh,Vh) = (vr^-u, iThv). Then, the first three terms on 
the right hand side can be estimated with the help of Lemma l4.1l and Holder's 
inequality: 

f(u-7T h u) < |H| wM( n) (E4 /3 / I/I 473 )' 74 , 

J2 [WE-<r h ](ti-ir h u)\ <M w l2 (m) ( E h E [\We ■ cr h ]\ 2 y 2 
Enan=<b JE Enan=0 Jh 

and 

\(Sh(uh\an + v h - no) + {"an ■ <?h) ~ *o)> (« - n h u)\ 9 n + v - n h v)\ 

< (\\u\\ l, — I — 1 1 "0 1 1 l, + \\Vh II 1, ) \\Sh(uh\fin + Vh — Un) + (van • <Jh) — tn\\ ., 

Note that the boundedness of the set of minimizers, Lemma 12.11 provides 
an explicit uniform bound on both \\u, v\\x an d \\uh, Vh\\x i n terms of the norms 
of the data. The W^ - 2 ,2 (<90)-norm leads to r/s ([SJ). 

The remaining term requires a slightly more precise analysis. Decompose 

("an ■ <?h) = {van • &h)+ - {van • Vh)- 

into its positive and negative parts. For a classical exact solution, the Signiorini 
condition requires (vqq ■ o) + = 0, and we estimate the corresponding term as 
above: 

1/2 
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For the negative part, we may drop the unknown term: 




The a posteriori estimate follows. □ 
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